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Abstract. We investigate the scalar Kir 
form factor at low energies by the method 
of unitarity bounds adapted so as to in- 
clude information on the phase and mod- 
ulus along the elastic region of the uni- 
tarity cut. Using at input the values of 
the form factor at t — and the Callan- 
Treiman point, we obtain stringent con- 
straints on the slope and curvature pa- 
rameters of the Taylor expansion at the 
origin. Also, we predict a quite narrow 
range for the higher order ChPT correc- 
tions at the second Callan-Treiman point. 

1 Introduction 

The low energy properties of the Kir form factors are of 
great interest both experimentally and theoretically. In 
particular, a precise knowledge of the slope and curvature 
parameters at t = would serve to improve the exper- 
imental analysis of K^ decays, confirm the predictions 
of chiral perturbation theory (ChPT) and provide bench- 
marks for future lattice determination of these quantities. 

In the present paper we consider the scalar Kir form 
factor, expanded as 
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in the physical region of K^ decay. The dimensionless pa- 
rameters Ag and Xq are related by X' Q = M%(r^ K )/6 and 
Aq = 2M*c to the radius (r^ K ) and curvature c used al- 
ternatively in the literature. 

The Kir form factors have been calculated at low ener- 
gies in ChPT [TJUMI an d on the lattice (for recent reviews 



see HE]). At t = 0, the present value /+(0) = 0.962±0.004 
[4] shows that the corrections to the Ademollo-Gatto theo- 
rem are quite small. Other low energy theorems frequently 
used are [5]-[l| 
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where Ak-k = M^ — M£ and Ak-k = —Ak-k are the first 
and second Callan-Treiman points, respectively. The low- 
est order values are known from Fk/F^ = 1.193 ± 0.006 
[I], and the corrections calculated to one loop are Act = 
-3.1 x 10~ 3 and A C t = 0.03 [2]. The higher order correc- 
tions appear to be negligible at the first point jS], but are 
expected to be quite large at the second one. 

Analyticity and unitarity represent a powerful tool for 
obtaining information on the Kir form factors. Several 
comprehensive dispersive analyses were performed recently, 
using either the coupled channels Muskhelishvili - Omnes 
equations [9"lll0j , or a single channel Omnes representation 

Alternatively, the method of unitarity bounds, pro- 
posed a long time ago in |121ll3j. and applied since then 
to various electromagnetic and weak form factors, exploits 
the fact that a bound on an integral of the modulus squared 
of the form factor along the unitarity cut is sometimes 
known from independent sources. Standard mathematical 
techniques then allow one to correlate the values of the 
form factor at different points or to control the truncation 
error of power expansions used in fitting the data [14| . 

For the Kir system the method was applied in [T3J and 
more recently in |16(I17| . In Ref. [TB] the method was ex- 
tended by including the phase of the scalar form factor 
along the elastic part of the cut, known from the elastic 
Kir scattering by Watson's theorem, while in [T7] infor- 
mation on the form factor at the second Calln-Trciman 
point was included for the first time in the frame of the 
standard bounds. 

In the present work we revisit the issue of bounds on 
the expansion coefficients ([!} by applying a more sophis- 
ticated version of the unitarity bounds proposed in |18) . 
The method uses the fact that the knowledge of the phase 
allows one to remove the elastic cut and define a func- 
tion with a larger analyticity domain. To be optimal, the 
method requires also some information on the modulus of 
the form factor in the elastic region. In [16], where the 
phase constraint was treated using Lagrange multipliers, 
this stronger property of the phase was not exploited, since 
no experimental information on the modulus was available 
at that time. 

More recently, the precise measurements of the r — > 
K-kv t spectral function by Belle collaboration [THj pro- 
vided also a first direct experimental determination of the 
modulus of the Kir form factors below a certain energy. 
The modulus is available also from the dispersive anal- 
yses [91HU], This justifies the application of the method 
proposed in |18j . Our work extends the analysis made in 
|17) by including information on the phase and modulus of 
the form factor on a part of the cut, which leads to a con- 
siderable improvement of the bounds. In the next section 
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we describe briefly the method for the scalar form factor, 
and in section [4] we present the results. A more detailed 
analysis, including a discussion of the experimental impli- 
cations and of the vector form factor, will be presented in 



2 Standard and new unitarity bounds 

The method makes use of the following mathematical re- 
sult: let g(z) be a function analytic in the unit disk \z\ < 1 
of the complex z-plane, which satisfies the inequality: 
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where / is a positive number. Then, if 

g(z) = .9o + giz + 92Z 2 + 
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is the Taylor expansion of g(z) at z = 0, and g(z±), g(z%) 
denote the values of g at two points inside the analyticity 
domain, \zi\ < 1, \z 2 \ < 1 (for simplicity we assume that 
Z\ and Z2 are real), the following determinantal inequality 
holds: 
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Moreover, all the principal minors of the above matrix 
should be nonncgative. For the proof see Refs. [1211131115] . 

To obtain this formulation for the scalar Ktt form fac- 
tor, one starts from a dispersion relation for the scalar 
polarization function of the s and u quarks [T5l[T6lT7ll2"l"] : 
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where unitarity implies the inequality: 
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with t± = (Mr ± M v ) 2 - We use here the notations from 
[21"] . where 7T is defined as the longitudinal part of the 
correlator of two vector currents. As in [IB], II o can be 
identified with (<?(f) — &(0))/t 2 , where is the correlator 
of the divergence of the vector current. 

The quantity x (Q 2 ) m @ can be reliably calculated 
by pQCD when Q 3> Aqcb- At present, calculations avail- 

~ [13] give: 
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where the running quark masses and the strong coupling 
a s are evaluated at the scale Q 2 in M S scheme. 



Taking into account the fact that fo(t) is analytic ev- 
erywhere in the complex i-plane except for the branch 
cut running from t+ to oo, the relations ([S])-© can be 
expressed in the canonical form ([3]) if one defines the vari- 
able 



t 



(9) 



which maps the t plane cut from t + to oo onto the unit 
disk \z\ < 1, such that z(0) = 0, and the function 

g(z) = f (t(z))w(z), (10) 

where t(z) is the inverse of © and w(z) is the outer func- 
tion [IS]- [HQ: 
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It may be noted that z is an independent variable in the 
outer function, whereas z(—Q 2 ) etc., are defined via the 
conformal variable eq.(j9]). 

We use now the fact that, below the inelastic thresh- 
old £i n , the phase of the form factor fo(t) is known from 
Watson's theorem and the 1 — 1/2 S-w&ve of elastic Ktt 
scattering. Then one can define the Omnes function 



0(t) = exp 



t 



6(f) 



dt'- ,, 



(13) 



where 5(t) is the phase of the form factor known for t < 
tin, and is an arbitrary function, sufficiently smooth (i.e. 
Lipschitz continuous) for t > t- m . It can be shown [20] that 
the results are independent of the function 5(t) for t > t- m . 

Since the Omnes function 0(f) fully accounts for the 
second Riemann sheet of the form factor, the function h(t), 
defined by 

f (t) = h(t)O(t), (14) 

is real analytic in the t- plane with a cut only for t > t- m . 
Then, the relations (0-© and (fT4]l can be expressed in 
the canonical form ©, by defining the new variable [T8] 
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which maps the t-plane cut for t > tj, 
\z\ < 1 of the z-plane, such that z(0) 
function HH] 



(15) 



onto the unit disk 
= 0, and define the 



g(z)=f (t(z)) W (z)u(z) P(t(z))} 



(16) 



x We mention that a \/2 is missing in the corresponding 
expressions given in |16II17| . 
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where t(z) is now the inverse of (|15[) . The new outer func- 
tion w{z) is defined as 
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and is calculable if the modulus |/o(£)| is known at low en- 
ergies, below £; n . Thus, we can use the inequality ([5]) and 
the nonnegativity of the leading minors to obtain bounds 
on the parameters of the expansion (JTJ. The Taylor coef- 
ficients in are defined uniquely in terms of these pa- 
rameters by (ITU|) or (|16|) . We further choose z\ = z(Ak-k) 
and z 2 = z(Aktt), where z is defined by ([9]) or (|T5|) . and 
express g(zj) in terms of the values in by using either 
(HHD or 

In our analysis we take as inputs the values of the 
form factor at £ = and t = Ak-k, the phase below i; n 
and the integral over the modulus required in (|19j) . Then 
the constraints resulting from ([5]) restrict the coefficients 
Aq , Aq and the value of fo (£) at the second Callan-Treiman 
point Ak-k- 



3 Input 

We work in the isospin limit, adopting the convention that 
Mk and M n are the masses of the charged mesons. The 
inputs provided by the low energy theorems was discussed 
in the Introduction. 

For choosing t; n , we recall that the first inelastic thresh- 
old for the scalar form factor is set by the Krj state, which 
suggests to take t- m = 1 GcV 2 as in [16]. However, this 
channel has a weak effect, the elastic region extending 
practically up to the K 7/ threshold, which justifies the 
choice ti n = (1.4GeV) 2 . In our analysis we shall use for 
illustration these two values of i; n . 

Below tin the function 8{i) entering (jl3|) is the phase 
of the S'-wave of I = 1/2 of the clastic Kir scattering 
[2"4"1ITU] . In our calculations we use as default the phase 
from [TO]. Above t- m we assume S(t) as a smooth function 
approaching ir at high energies. We checked numerically 
that the bounds are independent of the choice of <5(t) for 

£ > t in . 

To estimate the integral appearing in (|19[). we first 
used the parametrization of |/o(£)| in terms of the reso- 
nances k and -ftTo(1430), proposed by Belle collaboration 



j!9| . Using as input the solution 1 in Table 4 and Eq. 
(7) of [19], the integral has the value 66.08 x 1(T 6 for 
tin = lGeV 2 , and 184.89 x 1CT 6 for t in = (1.4 GcV) 2 . 
Although the parametrization used in j!9) does not have 
good analytic properties, this fact is not relevant for our 
analysis: all that we need is a numerical estimate of the 
integral in (fT9"|) . The analyticity of the form factor is im- 
plemented rigorously in our approach, for every numerical 
input. 

Aternativcly, using the modulus available from the dis- 
persive analyses [21] or [TU], the low energy integral in 
([T9|) is 40.05 x 10" 6 or 37.01 x 10" 6 , respectively, for 
tin = 1 GeV 2 , and 89.31 x 10~ 6 or 81.26 x 10~ 6 for t in = 
(1.4 GcV) 2 . 

Finally, we take Q 2 = 4 GeV 2 as in [TS P T ] . and ob- 
tain xo = (253 ± 68) x 10" 6 , using in ® m s (2GeV) = 
98±10McV,m u (2GcV) = 3 ± 1 MeV g] and a s (2GcV) = 
0.308±0.014, which results from the recent average a s (m T ) 
0.330 ± 0.014 (2M26][27l[28]. The error of xo includes also 
a contribution of 15%, of the order of magnitude of the 
last term in (O, to account for the truncation of the ex- 
pansiorQ. 



4 Results 

In order to illustrate the effect of the additional informa- 
tion on the phase and modulus, we compare in Fig. [T] the 
allowed domains in the plane (Aq, Aq), obtained with the 
standard and the new bounds, using only the constraint 
at t = (this case is obtained from ([5]) by removing the 
lines and columns that contain g(zi) and g(z2)). The large 
ellipse is obtained with the standard bounds, (j9])- (fl"2"j) . the 
small ones represent the new bounds, calculated with (|15[) - 
([T!?]) for two values of t in . For /o(0) we took the central 
value 0.962. The left panel is obtained with the integral in 
([T2"]) calculated with the modulus from [15] , for the right 
one we used the modulus from |24) . 

The inner ellipses are slightly smaller in the left panel 
than in the right one, because in the latter case the inte- 
gral in ([T9]) is smaller and / is larger (it is easy to see that 
a larger value of / leads to an ellipse of a larger size). 
Moreover, in the right panel the small ellipses are not 
contained entirely inside the large one, which means that 
among the functions satisfying the constraints (fl"5|) -(Jl"9" l) 
there are some that violate the original bounds ([9"j)- ([T2"|) . 
However, this does not mean that the conditions are in- 
consistent, since the ellipses have a nonzero intersection, 
which represents the allowed domain in this case. 

The effect of the low-energy theorems ^ is illustrated 
in Fig.[5J which shows the improved bounds (|15I) - ([T9")) cal- 
culated with the modulus from [TS] for ti n = 1 GcV 2 . The 
large ellipse is obtained using the constraints at £ = and 
t = Ak-k-i the tiny ellipses result from the constraints at 



2 In fact, it is known that the perturbative series in QCD 
are divergent. Improved expansions exploiting this feature were 
proposed, see for instance [2S]. However, in the present context 
this improvement is not necessary, since the sensitivity of the 
bounds to the value of xo is quite low. 
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X' x 1 



Fig. 1. Allowed regions for the slope and curvature from the 
standard and the new bounds with the constraint at t = 0. 
Left: modulus from [TS]; Right: modulus from [24] . 



t = 0, Ak-k and ^5^^, for the central values of /o(0) and 
Jo{Ak-k) and several values of the correction Act- The 
strong constraining power of the simultaneous constraints 
at Aktt and Akk was noted in |17| . However, the bounds 
derived now are much stronger than those in [17] . due to 
the additional information on the phase and modulus on 
the cut. 



A CT = -0.0134 



"* 0.6 - 



1.5 1.6 

X'x 1 2 



Fig. 2. Large ellipse: new bounds obtained with ti n = 1 GeV 2 
and the values at t = and Ak-k', small ellipses: new bounds 
for ti n = 1 GeV 2 and simultaneous constraints at t = 0, Ak-k 
and Ak-k, for several values of Act- 



The small ellipses exist only for Act inside a rather 
narrow interval, whose end points lead to inner ellipses of 
zero size in Fig.[3J Actually, this range results directly from 
the inequality by keeping only the lines and columns 
involving go, g(z{) and g{z%) and using the central values 
at t = and A k -k, we obtain -0.046 < A C t < 0.014. 

We recall that the current ChPT prediction is —0.057 < 
Act < 0.089 (cf. Eq. (4.9) of [TJ] adapted to our input 
Fk/Ftt)- As this interval is larger than the range derived 
above, we conclude that, at present, one can not further 
restrict the domain for the slope and curvature using the 
low-energy theorem at the second Callan-Treiman point: 
by varying Act inside its currently known range we ob- 
tain the union of the tiny ellipses in Fig. [31 which covers 
the large ellipse obtained using only the value at Ak-k- 

In Fig. [3] we show the allowed domains for the slope 
and curvature obtained with the constraints at t = and 
Ak-k for two values of i; n in the left panel. In the right 



panel we also superimpose the allowed region when no 
phase and modulus information is taken into account. As 
in Fig. [21 the low energy integral in (flT)| was calculated 
with the modulus from [19]. For t- m = 1 GeV 2 the large 
ellipse implies the range 0.0137 < X' Q < 0.0172, for i in = 
(1.4 GeV) 2 the small ellipse implies the narrower range 
0.0150 < Aq < 0.0163. In both cases we obtain a strong 
correlation between the slope and the curvature. It may 
be clearly seen that a dramatic improvement is obtained 
by the inclusion of phase and modulus data. 



— 1 m = 1.0 GeV 

— I. =(1.4) ! GeV 2 




w 



1.4 1.5 1.6 1.7 

X„x 10 2 



Fig. 3. Left: Allowed regions for the slope and curvature using 
as input the phase and modulus up to t\ n , and the values of 
/o(0) and /c^k,); right: as in left panel and also showing the 
region obtained with no phase and modulus information. 



As discussed, the method gives also very sharp pre- 
dictions for the corrections Act at the second Callan- 
Treiman point: for t- m = (1.4 GeV) 2 we obtain —0.031 < 
Act < -0.008 

The above ranges were obtained using the central val- 
ues of /o(0), /o(Ak-k) and xo, the phase from [TU] and the 
modulus from |19j . Accounting for the errors and using 
alternatively the phase and modulus from [24] . the end 
points of the range of Aq for t- m = (1.4 GeV) 2 varied by 
±0.00039| ± 0.00044^ ± 0.00028| Xo ± 0.00038| mod ± 
0.00070| ph , while for A C t the variation was ±0.0074| ± 
0.0002U^ ± 0.0052| XO ± 0.0069| mod ± 0.0020| ph . We note 
that, while the bounds are very sensitive to the input value 
of /o(0), the uncertainty of xo has a relatively low influ- 
ence on the results. 

In conclusion, our analysis shows that the modified 
type of unitarity bounds proposed in [Tg] , which includes 
input from the elastic part of the cut, leads to very strin- 
gent bounds on the scalar Ktt form factor at low energy. 
Using as input the precise values at t = and Ak-k and as- 
suming that the inelasticity is negligible below (1.4 GeV) 2 , 
we obtain for the slope at t = the range 0.0150(10) < 



X' < 0.0163(10), where the error is obtained by adding 
in quadrature the uncertainties due to various inputs. As 
shown in Fig. [31 the method leads to a strong correlation 
between the slope and the curvature. We obtain also a nar- 
row admissible range -0.031(12) < Act < -0.008(12) for 
the higher order ChPT corrections at the second Callan- 
Treiman point Ak-k, significantly reducing the range from 
ChPT mentioned earlier. Unlike in the usual dispersive 
approaches, the predictions arc independent of any as- 
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sumptions about the presence or absence of zeros, or the 
phase and modulus of the form factor above the inelastic 
threshold. 
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